Universal Distribution of Centers and Saddles in Two-Dimensional Turbulence 
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The statistical properties of the local topology of two-dimensional turbulence are investigated 
using an electromagnetically forced soap film. The local topology of the incompressible 2D flow 
is characterized by the Jacobian determinant A(x,y) — j(lo 2 — a 2 ), where u(x,y) is the local 
vorticity and a{x,y) is the local strain rate. For turbulent flows driven by different external force 
configurations, P(A) is found to be a universal function when rescaled using the turbulent intensity. 
A simple model that agrees with the measured functional form of P(A) is constructed using the 
assumption that the stream function, ip(x,y), is a Gaussian random field. 



Geometrical and topological ideas have played impor- 
tant roles in our understanding of turbulent phenom- 
ena. For example the energy cascades of both three- 
dimensional (3D) and two-dimensional (2D) turbulence 
are generally pictured as hierarchical structures of inter- 
acting eddies Q|. Although topological ideas find wide 
spread use in the interpretation of data and theory, they 
are not generally used to analyze data and develop the- 
ory. The reason for this is because there currently exist 
no measurement technique for the extraction of entire 
turbulent velocity fields from a 3D fluid. Since the bulk 
of turbulence experiments have been performed in 3D, 
this lack of enthusiasm for topological approaches is an 
artifact of limited data. In 2D, however, the situation is 
more tractable since particle tracking can measure entire 
2D velocity fields 0. In this article the statistical prop- 
erties of the local topology of 2D velocity fields extracted 
from a turbulent soap film are investigated. 

The local flow topology in a 2D fluid can be 
determined, up to a coordinate transformation, by 
two numbers: the Jacobian determinant, A(x, y) = 
d(v x ,v y )/d(x,y), and the divergence of the flow field, 
V - v. Since the velocity fields under consideration are es- 
sentially incompressible, the single field A(x, y) uniquely 
defines the flow. This quantity lends itself to a simple 
geometrical interpretation. If A(x, y) < 0, then in the 
reference frame of the fluid at point (x, y) the flow struc- 
ture is similar to that shown in the inset of Fig. 0(a) and 
is called a saddle. On the other hand, if A(x,y) > 0, 
the flow structure resembles that in the inset of Fig. l](b) 
and is called a center. If A(x,y) = 0, the local field is 
in either a state of pure shear, or more complicated non- 
linear shear. The quantity A also plays a significant role 
other than measuring the local topology of a 2D fluid. 
If the 2D fluid has density p and experiences external 
forces that are divergence free, then A is proportional to 
the Laplacian of the pressure field, Ap = 2pA. Therefore 
saddles are regions of locally high pressure, while centers 
are regions of locally low pressure. 

One can express A in terms of two more familiar quan- 



tities, the vorticity lo and the strain rate a: 

A= X -y-a\ (1) 

where uj 2 = Y,, ,['>,>' , - djVi) 2 /2 and a 2 = J2i,j( 9 i v j + 
djVi) 2 /2. Hence, regions of large positive A correspond 
to strong vorticity, whereas regions of large negative A 
correspond to strong elongation. Note that this is not 
exclusive, meaning that saddle regions could also have 
a large vorticity provided the local strain is still larger. 
Similarly, centers could be strongly elongated as long as 
their rotation rate exceeds the strain rate. 

This article will investigate the probability P(A) that 
a point in turbulent flow is behaving as a saddle or a 
center of a given strength. This will be done in 2D us- 
ing an electromagnetically forced soap film (e-m cell) [[| . 
Theoretical predictions of the behavior of P(A) are also 
developed in this article and shown to be in agreement 
with most of the features observed in the measurement. 

A description of the operation of the e-m cell can be 
found in |3|] and will not be reviewed in depth here. 
Briefly, the e-m cell creates a state of energetically steady 
turbulence in a 50 /an thick soap film using electromag- 
netic fields. Data from the e-m cell has been shown to be 
consistent with exact predictions that can be made using 
the forced 2D incompressible Navier-Stokes equation: 
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Here, v is the velocity field, p is the internal pressure 
field normalized by the fluid density, and F is the elec- 
tromagnetic force field acting on the fluid. The constant 
v ~ 0.016 cm 2 /s is the kinematic viscosity. The last 
term in Eq. ||, — av, accounts for the frictional effect of 
air acting on the soap film. This system has two di- 
mensionless control parameters which can be varied by 
changing either the magnitude of F or a. The first is 



the Reynolds number defined as Re 



s/ki n jV and 



the second is the dimensionless air- friction 7 = a/kf^v. 
Here ki n j is a typical wavenumber describing the spatial 
variation of the electromagnetic force field, called the in- 
jection wavenumber. In these experiments, 7 assumes 
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FIG. 1: A field in turbulence driven by a linear array of mag- 
nets. For clarity, A~ and A + are plotted in (a) and (b) re- 
spectively with darker shades indicating larger magnitudes of 
A* . The solid lines are the contours of the stream function, 
tp, and the large tick marks correspond to distances of 1 cm. 
The insets in (a) and (b) show the local flow for a saddle and 
center, respectively. 



values from 0.1 to 1, corresponding to a regime of weak 
to moderate damping, and Re varies from 10 to 100. 

To establish that the local topology measurements 
are not merely a reflection of the geometrical aspects 
of the external forcing two different force configurations 
were used. Both force fields were unidirectional with 
F y = 0. The F x for each field varied in different spa- 
tial patterns which will be called linear and square. F x 
is given by F x — Fq sin(fcoj/) for the linear field, and 
F x = Fq sin(fco2//v / 2) sin(fco^/\ // 2) for the square field. 
For both cases, ki n j = ko = 27r/a, with a ~ 0.6 cm for 
the linear forcing and a ~ 0.8 cm for the square forc- 
ing. Unless otherwise noted the magnitude of the force 
fields, Fq, varied in time with a square waveform of 3 Hz. 
Finally, topological aspects of decaying turbulence were 
also investigated by driving the cell and suddenly switch- 
ing off the drive. It is generally believed that decaying 
turbulence is very different from forced turbulence due 
to the absence of an inverse energy- cascade range jlj . 

The velocity field is measured by particle tracking ve- 
locimetry [|| , which typically yields 5000 velocity vectors 
per image. Figure |l| is a typical A field calculated from 
a single velocity field driven by the linear forcing with 
v rms ~ 10 cm/s. For clarity, the distribution of centers 
A > 0, denoted as A + , is plotted separately from the dis- 
tribution of saddles A < 0, denoted as A~. Unlike the 
velocity fields, which tend to be spatially extended, the 
distribution of A appears to be localized or spotty, more 
so for centers than for saddles. 

Sequences of 400-500 A fields were acquired for various 
values of 7 and Re, and for different types of forcing (see 
Table |). From these, probability distribution functions 
(PDF) of A were computed, and the results are shown in 
Fig. ||. All of the PDFs have a strikingly similar form. 
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t (s) 
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symbol 


linear 


60 


0.28 


na 


0.31 


712 


2.53 





linear 


60 


0.56 


na 


0.31 


1014 


2.60 


D> 


linear 


60 


0.97 


na 


0.31 


1282 


2.54 





square 


80 


0.6 


na 


0.4 


898 


2.84 


* 


decay 


70 


0.5 


0.5 


0.35 


104 


2.54 


< 


decay 


70 


0.5 


1.0 


0.35 


39.1 


3.4 


A 



TABLE I: Global constants for different types of forcing in the 
e-m cell. Here, r is the amount of time the system was allowed 
to decay before taking data. The symbols correspond to those 
plotted in the later figures. Re and 7 values listed for decaying 
turbulence denote the initial values of these parameters. 




FIG. 2: (a) The distribution function P(A) for: (a) linearly 
forcing with Re ~ 60 and 7 values of 0.28 (o), 0.56 (>) and 
0.97 (o), (b) square forcing with Re » 80 and 7 ~ 0.6, and 
(c) decaying turbulence at two separate times, 0.5 s (<) and 
1.0 s (A) after the forcing is turned off. All of these data sets 
are plotted together in (d) after normalizing by the respective 



One observes that the PDF is asymmetric with a nearly 
exponential decay for A < and an initially quick, per- 
haps stretched exponential, decay for small A > 0, even- 
tually relaxing to what appears to be exponential decay 
for large A > 0. The anomalously long tail of the PDF 
for A + is consistent with the visual appearance of veloc- 
ity fields in that powerful vortices are more conspicuous 
than saddles in the flow. Finally, we note that the PDF 
is non- analytic at A = 0. 

The PDFs in Fig. ^(a)-(c) can be collapsed on a single 
master curve if A is normalized by its RMS value, x = 
A/A rms . As seen in Fig. ||(d), over four decades in the 
vertical scale, the quality of the data collapse is excellent. 
The agreement between decaying and driven turbulence 
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also indicates that P(A) is independent of the presence 
or lack of an inverse cascade. These observations lead us 
to believe that the functional form of -P(A) is universal 
for 2D turbulence, though the flow appears to be visually 
different. 

The asymmetry of P(A) is clearly reflected in the 
moments of the PDF. All measurements show that 
(A) w 0, yielding (A~) w — (A + ). However, the sec- 
ond moment is tilted toward the positive side with ( = 
((A+) 2 )/((A-) 2 ) - 2.5 (see Table @). Though the data is 
not presented here, it should be noted that these results 
in the turbulent regime are markedly different from those 
obtained in the laminar flow regime with an ordered ar- 
ray of vortices. In this case, the PDF is nearly symmetric 
with (A~) = -<A+) and <(A+) 2 ) « ((A~) 2 ). 

We also considered the conditional probability of en- 
strophy u> 2 and square strain rate a 2 for a given A. We 
find that the conditional PDF for both cu 2 and a 2 are 
approximately exponential. The conditional expectation 
values for (a 2 ) and (cu 2 ) for a given A are plotted in Fig. 
||. It is noticeable that while (a 2 ) and (ui 2 ) depend lin- 
early on A for A < 0, such a linear relation does not hold 
for small positive A. 

The universal nature of P(A) indicates that the prob- 
ability distribution can be understood from basic prin- 
ciples and minimal models without requiring a detailed 
model of the flow. Since the experiment approximates 
2D incompressible flow, the velocity can be written in 
terms of a stream function ip(x, y) with v x — d y rp and 
v y = —d x ip. The Jacobian determinant is then A = 

-{d X dy^) 2 + (d 2 ^(d 2 ^). 

We can show that (A) = if the correlation 
function of ip is translationally invariant, that is, 
(ip(x, y) tjj(x' ' , y')) = M(u, w) is a function of u = \x — x'\ 
and w = \y — y'\ only. If we further assume that M(u, w) 



is analytic at u 

(A) = « 



w 



0, we find 



-{d x v x )(d y Vy) + (d x v y )(d y v x )) 
■(d u d w d w d u M) + {d 2 u d 2 w M) = 0. 



(3) 



Note that all derivatives of M here and below are eval- 
uated at u = w = 0. We expect that (A) = for a 
wide variety of flows since we do not assume that M is 
isotropic, simply that it is translationally invariant. 

To understand the PDF in the turbulent regime, we 
need a further assumption. The simplest assumption is 
that the stream function is a Gaussian random field. The 
derivatives of a Gaussian field are also Gaussian so that 
d x v x , dyVy, d v v x and d x v y are all Gaussian random fields 
with their cumulants given by 4th order derivatives of 
M(u,w). Hence it is straightforward to calculate -P(A), 
since A is a functional of Gaussian fields. 

The PDF, -P(A), is found to depend on only two 
positive parameters. The first parameter, A = 
V((dxvy) 2 )((dyv x ) 2 ) = y/(d£M) tftM) simply sets the 
scale of A and hence does not affect the shape of the 
distribution. The second parameter 6 is 



b = - 



((d X Vy)(dyV X )} 



d 2 d 2 M 



V((d x v y ) 2 ) ((dyV x ) 2 ) y/(dtM) (0*M) 



(4) 



Hence b is the normalized correlation of the shear strain 
rates in the different directions and is the only parameter 
that affects the shape of P(A). In principle b is a free 
parameter. However, if M(u,w) has elliptical symmetry, 
that is M(u,w) = M(u 2 + cw 2 ) where c is a constant, b 
assumes the value 1/3. 

Introducing a rescaled A via A' = A/ A, the PDF can 
be written in terms of an integral Q : 



J 



P(A') = Cexp 



-A' 
l + b 



dr 



VA' + r 2 

r 



exp 



1-6 



2b 



(5) 



where r min — if A' > and r min = \/ — A' if A' < 0, 
and C is a normalization constant. For general b, this 
integral must be done numerically. However, for 6=1/3 
we obtain a close form: 

( 2Ctt 3 /2 

P(A') = 



3A'/2 



for A' < 0, 



2CV^ e 3A'/2 $ 



A' for A' > 0, 



(6) 



where $ c is the complemetary error function. The dif- 
ferent lower limits of integration, r m i n , for positive and 
negative A' makes P(A) asymmetric and non-analytic at 
zero. For negative A', P(A') is exponential while for pos- 
itive A', P(A') has a non-exponential behavior for small 
A' followed by an exponential decay at large A'. The 
asymptotic exponential decay rates k± can be obtained 



for general 6 with k + = 1/(1 + 6) and fc_ = 1/(26). Since 
6 < 1, k_ > k + and the asymptotic decay is faster for 
negative A'. 

To test this model, 6 was measured directly from the 
experimental flow fields and found to be close to the ex- 
pected value of 1/3; ranging between a low value of 0.3 
for linear forcing to a high value of 0.4 for the square 
forcing (see Table |). The scaled PDF determined from 
Eq. H for 6 = 1/3 together with the experimental PDFs 
are displayed in Fig. ||. The assumption that the stream 
function is a Gaussian field captures all the main features 
of the experimental PDF. In particular, the theory repro- 
duces the asymmetry between positive and negative A, 
and the exponential form of P(A) for negative A. For 
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FIG. 3: The conditional expectation values of a 2 (o) and J 2 
(•) for a given A. The measurement is taken from a linearly 
forced flow at Re = 60 and 7 = 0.56. 




FIG. 4: The same rescaled PDF as in Fifr 2(d) together with 
the theoretical PDF calculated from Eq. rawith b = 1/3 (solid 
line) . 



A > 0, the theory predicts the pronounced curvature on 
the semi- log plot in agreement with the experiment. For 
larger A the theory predicts that the PDf decays more 
slowly for positive than for negative A in agreement with 
observations. However, the experiment shows a much 
more pronounced curvature in the large positive A tail 
than the Gaussian model. This may indicate a break- 
down of the Gaussian assumption in regions of large vor- 
ticity. 

To summarize, the PDF of the Jacobian determinant in 
2D turbulence deviates strongly from Gaussian behavior 
with prominent exponential tails for large |A|. Moreover, 
P(A) is non-analytic at A = and decays more slowly 
for positive A or the centers. This is evidenced by the 
ratio ((A+) 2 )/((A~) 2 ) ~ 2.5. The non-analytic nature 
stems from the fact that centers and saddles are topolog- 
ically distinct and cannot transform continuously from 
one to the other. The experimental findings reported 
here are universal, independent of turbulent intensity and 
the means of turbulence generation. That this behavior 
is universal may be linked to the fact that A is a lo- 
cal quantity and therefore insensitive to the long-range 
spatial correlations of the turbulent velocity field. The 
locality of A may also explain the success of the Gaus- 
sian approximation for the stream function ip(x 7 y). This 
Gaussian assumption may prove useful in investigating 
other features of 2D turbulence including 2D scalar tur- 
bulence. 
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